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Magnetic excitations in the quasi-one-dimensional quantum antiferromagnets PbNi2V208 and 
SrNi2V2 08 are measured all the way up to the zone boundary energy using inelastic neutron scatter- 
ing from powder samples. An estimate for next-nearest-neighbor in-chain interactions is obtained. 
The role played by these interactions in spin-vacancy induced magnetic ordering in PbNi2V20g is 
discussed. 
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I. INTRODUCTION 

The quasi-one-dimensional (quasi-lD) S = 1 anti- 
ferromagnet (AF) PbNi2V208 was recently skown to 
be the first example of a Haldane-gap systemfil to un- 
dergo long-range magnetic-, ordering upon doping with 
non-magnetic impurities. ErH This unique behavior was 
understood through inelastic neutron scattering mea- 
surements, that suggested that three dimensional (3D) 
inter-chain interactions in PbNi2V208 are almost (but 
not quite) strong enough to destroy the HalcLaiie singlet 
ground state even in the undoped material.aa In fact, 
in the very similar isostructural compound SrNi2V208, 
inter-chain coupling exceeds the critical value, and, un- 
like the Pb-based material, the liBjdoped Sr-Ni vanadate 
orders magnetically at T = 7 K.Ha 

The introduction of non-magnetic vacancies into 
gapped quantum spin chains produces free S = 1 /^-de- 
grees of freedom on either side of each impurity sitcBij In 
the presence of arbitrarily weak inter-chain interactions, 
these free spins will order in three dimensions at suffi- 
ciently low temperature. However, next-nearest neighbor 
(nnn) intera|Ctions within the chains can prevent this from 
happening Naively, in the presence of AF nnn coupling, 
one can expect the liberated spins around each impurity 
to bind into non-magnetic singlets (dimers), preserving 
the disordered ground state of the system. 

There are at least two indications to that nnn inter- 
actions may be active in PbNi2V20g and SrNi2V20g. 
First, the spin chains in both materials are not linear, 
but are formed by Ni^"*" ions arranged in 4-step spirals 
propagating along the c axis of the orthorhombic lat- 
tice (for a detailed discussion of chain geometry and 
inter-chain interactions see Ref. In this "twisted" 
structure nnn superexchange pathways involving equa- 
torial oxygens of the NiOg tetrahedra can be quite effec- 
tive. Second, there remains some discrepancy between 



the in-chain exchange constant J deduced from the Hal- 
dane gaps measured with neutron scattering, assuming 
only nearest-neighbor (nn) interactions (J — 9.0 meV), 
and from high-temperature {T w J) susceptibility data 
(J = 8.2 meV). The former technique probes low-energy 
excitations around g|| = tt (AF zone-center), while the 
latter is sensitive to excitations in the entire Brillouin 
zone. The mismatch may be interpreted in terms of 
an additional coupling with a characteristic interaction 
length of twice the spin-spin separation in the chains, 
that is known to have different effects on zone-center and 
zone-boundary excitation energies. 

The most direct way to resolve this issue would be 
to measure the dispersion of magnetic excitations in the 
entire Brillouin zone in a single crystal sample. Unfor- 
tunately, no single-crystal or aligned powder samples of 
either PbNi2V208 or SrNi2V208, suitable for inelastic 
neutron scattering studies, are available to date. As will 
be discussed below, some information on nnn coupling 
can be obtained using powder samples, if the measure- 
ments are extended all the way to the zone-boundary 
energy. Previous neutron scattering studies employed 
cold neutrons to inwjatigate low-energy excitations in 
the two compounds In the present paper we report 
thermal-neutron studies specifically aimed at measur- 
ing high-energy excitations. We then engage in a semi- 
quantitative discussion of the role nnn interactions play 
in vacancy-induced ordering of weakly-coupled Haldane 
spin chains, and apply the results to the particular case 

of PbNi2V208. 



II. EXPERIMENTAL PROCEDURES AND 
RESULTS 

Both vanadates crystallize in tetragonal structure, 
space group /41cd, with lattice constants a — 
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12.249(3) A, c = 8.354(2) A for PbNiaVaps,^ and 
a = 12.1617 A, c = 8.1617 A for SrNisVzOgjB respec- 
tively. All measurements were performed on large (« 5 g) 
powder samples of PbNi2V208 and SrNi2V20|a, similar 
to those used in previous cold- neutron studies The ex- 
periments were carried out at the HB3 3-axis spectrom- 
eter installed at the High Flux Isotope Reactor at Oak 
Ridge National Laboratory. Neutrons with a fixed final 
energy of 13.5 meV were used with a 5 cm thick Py- 
rolitic Graphite (PG) filter positioned after the sample 
to eliminate higher-order beam contamination. PG(002) 
reflections were used for monochromator and analyzer, in 
combination with 48' - 40' - 40' - 120' coUimators. Sam- 
ple environment was in all cases a standard He-4-flow 
cryostat. The data were taken at T = 1.5 K. 
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FIG. 1. Typical background-corrected constant-Q scans 
measured in PbNi2V208 at T = 1.5 K. Arrows roughly indi- 
cate the gap and zone-boundary energies. Solid lines are sim- 
ulations based on the result of a semi-global fir to the data, 
that allows for next-nearest-neighbor in-chain interactions, as 
described in the text. Dashed lines are a best fit assuming 
only nearest-neighbor interactions are present. Statistical er- 
ror bars are smaller than symbol size. 

Typical data sets collected in constant-Q scans for 
PbNi2V2 08 and SrNi2V2 08 are shown in Fig. |] and 
Fig. 1^, respectively. Figure |3| shows the bulk of the data 
collected for the Sr system (only 5 scans, at |Q|/c* = 
2.5, 3, 3.5, 4 and 5, were collected for the Pb-based ma- 
terial). The measured intensities were scaled to correct 
for the previously determined energy-dependent A/2 con- 
tribution to the incident beam, that affects the monitor 
rate. Within the incident neutron energy range used. 



the monitor efficiency changes by roughly 30%. In ad- 
dition, a flat background (17 counts/min), measured at 
30 meV energy transfer (above the magnetic excitation 
band), was subtracted from the data. 
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FIG. 2. Typical background-corrected constant-Q scans 
measured in SrNi2V208 at T = 1.5 K. Arrows and solid lines 
are as in Fig. |l| Statistical error bars are smaller than symbol 
size. 

For both materials all scans have the characteristic M- 
shape. In a powder sample one measures the spherical 
average of the dynamic structure factor S{Q,uj). For 
large momentum transfers, a large portion of Q space is 
sampled simultaneously, and the powder cross section is 
in many ways similar to the density of states function 
n(Lu). In particular, one expects to see a peak in the 
powder cross section wherever the density of states is 
large, i.e., at the gap eoergy and at the zone boundary. 
As discussed previously,Q the peak around 5 meV energy 
transfer corresponds to excitations close to the Haldane 
gap at a momentum transfer g|| — n + 2?™ (n-integer) 
along the chain axis. The weaker higher-energy peaks, 
represent excitation near the top of the band, i.e., in the 
vicinity of = n/2 + nn (n-integer). 

III. DATA ANALYSIS 

A. Model cross section 

The convoluted nature of inelastic powder data make 
it notoriously difficult to interpret (for an example see 
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Rcf. particularly Section IV). The only way to extract 
information on the dynamic structure factor from such 
measurements is to assume a "reasonable" model and 
refine the parameter values to best-fit the data. The 
main assumption of the model used for PbNi2V208 and 
SrNi2V2 08 to date is that these compounds are described 
as weakly-interacting S — 1 AF chains with magnetic 
anisotropy. This picture is based on susceptibility, high 
field and specific heat measurements, and is supported bjz. 
the observed unique shape of constant- neutron scans.u 
In the present work we shall adopt this model as well. 

The form of the neutron scattering cross section for this 
model was derived in Ref. ^, and is written in the sin- 
gle mode approximation (SMA). Inter-chain interactions 
were accounted for within the Random Phase Approx- 
imation (RPA). This cross section function performed 
rather well in describing low- energy neutron data. Pre- 
vious stoidies of other Haldane-gap compounds, such as 
NENPB and CsNiClgy, have demonstrated that the 
SMA is, in fact, applicable in most of the Brillouin zone. 
The previously derived SMA cross section (Eqs. 5-13 
of Ref. ^ could thus be safely used in interpreting our 
high-energy results. However, the expression for the dis- 
persion relation (Eqs. 14 and 15 of Ref. ^ had to be 
slightly modified. These modifications are discussed in 
the following paragraphs. 



1. Isolated chains with nearest-neighbor interactions 

In Ref. I the dispersion of excitations in non- 
interacting Haldane spin chains was assumed to be si- 
nusoidal and symmetric around the classical AF zone- 
boundary (7|| = 7r/2. In fact, since translational symme- 
try is not broken in the Haldane ground state, the ex- 
citation energy pjtrapolates to a larger value at 
than at g|| — *■ Tr.lla This effect can be taken into account 
empiricaliv by including an extra term in the dispersion 
relation.EJ The dynamic susceptibility of a single Haldane 
spin chain is then written as: 
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(1) 



(^..^ff = A|, + .2sin%||+a^cos^^. (2) 



In these formulas, which replace Eq. 22 in Ref. || App 
is the gap energy for the Haldane gap mode of a partic- 
ular polarization /?/3, is the static staggered suscep- 
tibility, and V is the spin wave velocity. The spin chains 
PbNi2V208 and SrNi2V208 contain four equivalent spins 
per crystallographic period along the c axis. Throughout 
this work the actual wave vector transfer will be denoted 
as Q = ha* + kb* + Ic* . The symbol q\\ = Qc/A in Eq. | 
will stand for the reduced wave vector transfer along the 
chains. 

The energetics of a Haldane-gap AF is by now very 
well established through numerical work by several 



authors. For an S* = 1 Heisenberg chain with single- 
ion anisotropy of type D{S^)^ the approximate values of 
parameters are as follows: 
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(A) « 0.41J, 
V « 2.49J, 
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(5) 
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The parameter a in Eq. g characterizes the asymmetry 
of the dispersion relation. Existing experimental data for 
NENPEland CsNiClay gives a = I Ah J « 0.58i; and a = 
I.IJ « 0.44w, respectively. The latter results inspires 
more confidence, as it was obtained with a much higher 
experimental resolution and gives a dispersion relation 
very similar to that |Obtained by direct diagonalization in 
finite- length chains .Ej Throughout the rest of the paper 
we shall assume a fixed at a = OAAv. 



2. Next-nearest-netghbor interactions 

Next-nearest-neighbor in-chain interactions, that we 
believe can be relevant for PbNi2V2 08 and SrNi2V2 08, 
can be accounted for within the RPA. Isotropic nnn cou- 
pling of magnitude J' gives an additional term in the 
r.h.s. of Eq. || that can be written as: 



^/3/3a2 j. 
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The primary effect of nnn interactions is to modify the 
gap energies: 
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(9) 



where A/3^.0 are the gaps in the absence of nnn cou- 
pling. The relative effect of nn interactions on the zone- 
boundary energy and spin wave velocity is considerably 
smaller, and can be ignored. To a good approximation, 
in analyzing the experimental data one can use the dis- 
persion relation ^, but treat the gap energies and v as 
independent parameters, rather than assuming the rela- 
tion H valid only in the case J' = 0. The ratio ^ = v/{A) 
can then be determined and used to obtain an estimate 
for J'. From Eqs. for J' < J one gets: 



C«eo(l-9.33y) 



(10) 



Here — vj (Aq) ~ 6.07 is equal to the equal-time cor- 
relation length in the absence of nnn coupling. 
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FIG. 3. The bulk of constant-Q scans measured in SrNi2V208 at T = 1.5 K. Solid lines are as in Fig. Statistical error 
bars are smaller than symbol size. 



3. Inter-chain interactions 

For the coupling geometry found in PbNi2V208 and 
SrNi2V2 08, switching on inter-chain interactions within 
the RPA gives an additional term 

- x^/^lpJi cos(^V2) [cos(^h) + cos(7rfc)] (1~^°'^^V2) 

(11) 

in the r.h.s. of Eq. ^. In this formula, which corresponds 
to the last terms in Eqs. 14 and 15 of Ref. ^, Ji is the 
magnitude of nearest-neighbor inter-chain interactions. 
For our present purposes, the possible small anisotropy 
of Ji will be ignored. 

B. Parameters for PbNi2V208 and SrNi2V208 

Our understanding of the physics of PbNi2V2 08 and 
SrNi2V2 08 is based on the knowledge of the follow- 
ing characteristics: the intrinsic gap energies A^^ and 
l^xx = Ayj, , the magnitude of inter-chain interactions Ji 
and the spin wave velocity v. Can these four parameters 
be unambiguously determined from the available experi- 
mental data? 

First, consider those quantities that characterize low- 
energy excitations near q\\ = tt. For PbNi2V208 the 
actual minima of the 3D dispersion, -Emin./J/S were in- 
dependently measured in high-jiield magnetization stud- 
ies of aligned powder sampleso'Q These energies are di- 
rectly related to A/3/3 and Ji. The third experimen- 
tal value need to unambiguously "untangle" these three 



quantities was probed mi previous high-resolution cold 
neutron measurements.Q Indeed, the well-defined lower- 
energy 5 meV peak in constant-g powder scans is located 
approximately at the top of the dispersion band perpen- 
dicular to the chains at = tt, and its position is di- 
rectly related to A/3/3 and Ji. From the combined high- 
field and neutron measurements A/3/3 and Ji can thus be 
reliably determined within the established model. For 
PbNi2V208 these parameters are: A^^ = 3.1(3) rijneV, 
A^:r = ^yy = 4.0(3) meV and Ji = -0.17(2) meVE 

For SrNi2V208, which is ordered at low temperature, 
the model of weakly-coupled Haldane spin chains is, 
strictly speaking, not applicable. Nevertheless, as the 
ordered moment is expected to be very small, the model 
is still a good approximation, except in the direct prox- 
imity of the magnetic Bragg peaks. From the analy- 
sis of the low-energy powder spectrum for SrNi2V208 
the following parameters were previously determined: 
A,^ = 2.8(4) meV A^^ = Ayy = 3.9(3) meV and 
Ji = -0.17(3) meV.ta 

TABLE L Parameters of the model Hamiltonian used to 
analyze the experimental data. 



Parameter PbNi2V208 SrNi2V208 

All 3.1 (fixed) 2.8 (fixed) 

A^^ = Ayy 4.0 (fixed) 3.9 (fixed) 

Ji -0.17 (fixed) -0.17 (fixed) 

a 1.1 (fixed) 1.1 (fixed) 

V 20.6(7) 21.2 (4) 
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The exact value of the remaining "high-energy" param- 
eter V has httle impact on the low-energy powder spec- 
trum, and was not treated as an independent variable in 
Ref. y. In contrast, v determines the high-temperature 
magnetic susceptibility of the system, that, provided 
fceT A, are insensitive to the details of the low-energy 
spectrum. This type of experiments yields v ~ 20.4 meV 
for PbNi2V208. The present study provides the most 
direct measure of the zone-boundary energy, and hence 
the spin wave velocity v, defined by the location of the 
higher-energy peak and cutoff in the constant-g powder 
scans. In analyzing these data we shall use the cross sec- 
tion defined in Eqs. 5-13 of Ref. |4| in combination with 
the dispersion relation given by Eqs. ^ and [l| above, us- 
ing only V as an adjustable parameter of the model. 

C. Fits to experimental data. 

In the data analysis, the powder- average of the param- 
eterized cross section was calculated numerically using a 
Monte Carlo routine. The result was further numerically 
convoluted with the Q — E resolution function of the 
3-axis spectrorfteter, calculated in the Cooper-Nathans 
approximation.113 Performing a true global fit to the data, 
using a single set of parameters to describe all measured 
scans, was not possible because of the scattering-angle 
dependence of neutron transmission in the sample. The 
effect is due to a large sample size and appreciable in- 
coherent cross sections of Ni and V nuclei. Transmis- 
sion in the sample was studied by measuring the inten- 
sity of incoherent elastic scattering, and found to change 
by « 30%, depending on scattering geometry. Within 
each constant- Q scan the correction is small and could 
be ignored. However, in analyzing the entire data sets, 
collected in a broad range of momentum transfers (scat- 
tering angles), a separate intensity prefactor had to be 
used for each const-Q scan. Thus, for each compound, 
the data were analyzed in a semi-global fit: a single ad- 
justable parameter v was used to simultaneously describe 
all measured scans, each with its own adjustable prefac- 
tor. The parameters were refined to best-fit the data 
using a least-squares algorithm. 

Excellent fits were obtained with v — 20.6(7) meV for 
PbNizVzOg and t; = 21.2(4) meV for SrNizVsOg. Simu- 
lations based on the refined parameter values are shown 
in solid lines in Figs. The fits correspond to ~ A.7 
for PbNiaVaOg and = 3.8 for SrNiaVaOg. Given that 
the overall shape of the measured scans is reproduced by 
our model quite well, the large values of are most likely 
due to the fact that experimental statistical error bars 
are less significant than unavoidable systematic errors, 
such as multi-phonon scattering in the sample. While 
such effects can lead to an inaccurate background sub- 
traction, they are not likely to infiuence the measured 
zone-boundary energy, as the latter is accurately pin- 
pointed by the well-defined higher-energy cutoff in the 



energy scans. The relevant parameters used to fit the 
data, both fixed and adjustable, are summarized in Ta- 
ble ?? for both systems. The refined values for v are to be 
compared to previous estimates based on the measured 
gap energies and Eq.^, assuming only nearest-neighbor 
interactions, v — 22.5(1.0) meV and v — 21.5|^.9) meV, 
for PbNi2V208 and SrNi2V208, respectively.a For com- 
parison, the dashed lines in Figs. ^ are simulations based 
on these latter value for PbNi2V208. 

It was important to verify that any uncertainty that 
may exist for the dispersion asymmetry parameter a does 
not have a significant impact on our determination of v. 
This was done by analyzing the PbNi2V2 08 data using a 
larger value a — O.SSi;, as previously experimentally de- 
termined for NENPJIj rather than a — 0.44u. Essentially 
the same quality of fits was obtained, and the refined spin 
wave velocity v — 19.9(8) meV is within the error bar of 
that obtained with a = 0.44u, as described above. This 
demonstrates that the data analysis procedure employed 
is fairly insensitive to dispersion asymmetry. 

We see that any deviations from the model involv- 
ing only nearest-neighbor interactions are small in either 
material. While for SrNi2V208, within the error bars, 
the discrepancy is undetectable, for PbNi2V208 it ap- 
pears to be more substantial. It is gratifying to see that 
for the latter system, the bandwidth measured in this 
work is in excellent agreement with estimates based on 
the high-temperature susceptibility (see above). As ex- 
plained in the previous section, the reduction of zone- 
boundary energy is a sign of antiferromagnetic nnn in- 
chain interactions. Using Eq. |l^ for PbNi2V208 we ob- 
tain J' w 0.4(2) meV. Detecting such small values of J' 
is only possible due to the large staggered susceptibil- 
ity of a Haldane spin chain (the gap is relatively small 
compared to v), and hence the large prefactor in Eq. |l^. 

IV. DISCUSSION 

To summarize, the previously proposed nearest- 
neighbor model for PbNi2V2 08 and SrNi2V2 08 should 
be an adequate description of these compounds in most 
cases, in the entire Brillouin zone. For PbNi2V208 the 
data suggest, barely outside the error bars, a small anti- 
ferromagnetic nnn in-chain coupling of a few tenths of a 
meV. But how strong should this coupling be to have an 
impact on long-range ordering in a doped material? We 
shall now discuss the general significance of nnn inter- 
actions to spin- vacancy-induced order in weakly-coupled 
integer-spin quantum chains, and then consider the par- 
ticular case of PbNi2V2 08. 

A. Interactions between end-chain spins 

The problem of spin-vacancy-induced ordering in cou- 
pled integer-spin quantum chains with only nearest- 



5 



neighbor in-chain interactions has been studied theoret- 
ically by Shender and Kivelson.Q The introduction of a 
non-magnetic impurity in place of a S* = 1 ion produces 
two free S ^1/2 spins located on the two ends of the sev- 
ered chain.ETEl In each of the two fragments the end-chain 
spin is delocalized on a length scale defined by the intrin- 
sic equal-time correlation length ^ « 6 of a Haldanc spin 
chain. In principle, existing 3D interactions in the sys- 
tem should lead to long-range magnetic ordering of these 
liberated end-chain spins at a non-zero temperature. 



(c) 
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FIG. 4. A schematic representation of magnetic inter- 
actions around a non-magnetic impurity site. End-chain 
5=1/2 degrees of freedom (a and d) appear to either side of 
the impurity site (b). For simplicity, the delocalized nature 
of these end-chain spins is not shown. The interaction Ji of 
a liberated spin (a) with a neighboring chain locally polarizes 
the latter, inducing a "pocket" of staggered magnetization 
(c). The induced magnetization, in turn, couples to the other 
free spin (d). The resulting effective coupling competes with 
direct next- nearest-neighbor in-chain interactions J'. 

The most interesting result of Ref. ^ is that if inter- 
chain interactions are below the threshold value at which 
even the undoped system orders in three dimensions, 
there must be a critical spin-vacancy concentration Xc 
below which the system remains disordered at T = 0. 
The mechanism of this behavior is illustrated in Fig. ^, 
which show spins in the vicinity of one non-magnetic im- 
purity. For simplicity, the delocalized nature of end-chain 
spins is not shown. By virtue of inter-chain coupling Ji 
the end-chain spin on one of the fragments of the severed 
chain induces a "pocket" of staggered magnetization on 
a neighboring chain. The size of the pocket is determined 
by the correlation length ^. The induced staggered mo- 
ment interacts (again via Ji) with the end-chain spin on 
the other fragment. This adds up to an effective 2nd- 
order coupling between end-chain spins of the order of 
jJ^ffl « zx-nJi, where z is the inter-chain coordination 
number. In Ref. || it is argued that this coupling is anti- 
ferromagnetic and will bind the two end-chain spins into 
a non-magnetic singlet. As a result, the system will re- 
main in a singlet ground state despite the introduction of 



a vacancy. Another way of looking at this phenomenon 
is that the effective nnn interaction J^g "patches" the 
broken Haldane chain at the impurity site with an anti- 
ferromagnetic bond, albeit a weak one. 

We believe that this argument is, in general, incorrect. 
Indeed, from Fig. ^ it is clear that, no matter what the 
sign of Ji is, J^ff is a,ctn&\\y ferromagnetic (J^jj < 0) and, 
if anything, binds the two end-chain spins into a mag- 
netic triplet, delocalized over the length scale ^. Thus, 
in the absence of nnn interactions, long-range ordering 
should occur for arbitrary small impurity concentration 
X. The ordering temperature will, of course, be exponen- 
tially small with xu 

This conclusion has to be modified in the presence of 
actual in-chain nnn AF interactions. For J' > \J'eff \ the 
arguments of Ref. |^ can be repeated verbatim to show 
that the system will remain disordered at T = for suf- 
ficiently small X. On the other hand, the doped system 
will order at a non-zero temperature for J' ^ l"^e//l- 
any given small concentration there will be a quantum 
phase transition from a disordered to ordered state at a 
critical strength of J' given by 

J'c ~ zx.Jl (12) 

These estimates are semi-quantitative at best, yet they 
provide a simple physical picture and emphasize the rel- 
evant energy scales. 

B. Application to PbNi2V208 

In PbNi2V208, that remains in a gapped state at 
T — !■ 0, long-range ordering can be induced by substi- 
tuting some, of the S = 1 Ni^+ ions with non-magnetic 
Mg^+ ions.l3 The transition was observed fot, all Mg- 
concentrations x ranging from 0.1 down to O.Ol.B Though 
the Neel temperature Tn was found to decrease rapidly 
with decreasing x, there is, to date, no clear evidence of a 
critical concentration x^ below which the system remains 
disordered. The fact that in PbNi2V208 doping-induced 
ordering is related to end-chain spins has been elegantly 
demonstrated by comparing the effects of Mg^+ (5 = 0) 
and Cu^"*" {S = 1/2) substitution on the Ni-sites. The 
latter is exactly half as effective, in terms of the effect on 
magnetic susceptibility, since the Cu-spin and the two 
end-chain spins become bound in a doublet state. In 
contrast with S' = impurity substitution, the introduc- 
tion of a S' = 1/2 impurity produces only one, not two, 
free spins in the systemJa 

To estimate the critical strength of nnn AF interac- 
tions that would interfere with long-range ordering in 
weakly-doped PbNi2V208, we have to estimate J^g. Us- 
ing the experimental values of J, and assuming z — 
2,1 from Eqs. |l| and |-@ we obtain = |J^fj| = 
—0.18(6) meV. Given the experimental error bars and 
the semi-quantitative nature of Eq. |l^, we can only con- 
clude that in PbNi2V208 J' is of the same order of mag- 
nitude as the critical value J^. The question of whether 



6 



there is a critical doping level below which PbNi2V2 08 
remains disordered, stays open. Hopefully, extending the 
Ty!{x) measurements to lower concentrations, more ac- 
curate single-crystal neutron scattering measurements of 
the gap energies, and the development of a quantitative 
theoretical model will resolve this interesting issue in the 
future. 



V. CONCLUSIONS 

i) This work provides a direct and accurate measure- 
ment of the zone-boundary energy in PbNi2V2 08 and 
SrNi2V2 08. The previously proposed model, based on 
weakly-interacting S — 1 quantum spin chains with 
nearest-neighbor interactions, is shown to be an adequate 
description of these compounds in the entire dynamic 
range of magnetic excitations, ii) Unlike thought previ- 
ously, weakly-coupled Haldane spin chains should order 
in three dimensions at arbitrary small concentrations of 
non-magnetic impurities, iii) Next-nearest neighbor AF 
interactions can inhibit long-range ordering at small im- 
purity concentrations, but only if they exceed a certain 
critical value J'^. iv) In the particular case of PbNi2V208 
nnn interactions, if at all present, are no larger that sev- 
eral tenths of a meV, and may be of the same order of 
magnitude as J^. 
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